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Abstract — We present a general method for the error analysis 
of spatial modulation (SM) systems over correlated and uncorre- 
cted Rayleigh and Rician fading channels. The proposed method, 
making use of the properties of proper complex random variables 
and vectors, provides an exact upper bound for the class of fading 
channels considered for any number of transmit and receive 
antennas and for a wide family of linear modulation alphabets. 
Theoretical derivations are validated via simulation results. 

I. Introduction 

Spatial modulation (SM) has emerged as a new and highly 
effective multiple-input multiple-output (MIMO) communica- 
tion technology Q), J2). The fundamental principle of SM 
is simple: Both the signal and the antenna spaces are used 
to convey information. In other words, at each transmission 
instant the incoming information bits are split into two groups. 
The first group is mapped onto a transmit antenna index 
whereas the second group is used to select a constellation 
symbol within the modulation alphabet. Then only the antenna 
designated by the antenna bits is activated to send the selected 
symbol. This approach has a number of advantages: i) Because 
one antenna is activated at each time instant, not only the 
inter-channel interference (ICI) is avoided but also the receiver 
complexity is reduced, ii) Given a spectral efficiency of R 
b/s/Hz., because SM allocates some bits for the antenna index 
instead of mapping all R bits to the constellation symbols, 
the constellation size is reduced in comparison to single 
antenna transmission. This increases the minimum distance 
between constellation points and results in an improved error 
performance, especially when R is large. Notice that, space- 
shift keying (SSK) proposed in [3| as a special case of 
SM where information is transmitted only over the antenna 
space, (i.e., SM with a constellation size of 1) have similar 
performance with even lower complexity, especially at low 
spectral efficiency. On the other hand, SSK is required to 
double the number of transmit antennas for every increment 
in the spectral efficiency. As a result SM is still a more viable 
alternative for practical number of antennas. 

Despite its advantages, the unique nature of SM systems, 
that is, the transmission and reception of two different forms 
of information simultaneously, creates a challenge in their 
understanding and analysis. One challenging aspect is the 
performance analysis of SM under optimal detection. Orig- 



inally a performance analysis is presented for Rayleigh fading 
channels in (TJ, which is shown to be suboptimal in 10. 
Then in [2], a maximum likelihood (ML) detection rule is 
presented along with a performance analysis for Rayleigh 
fading channels using the well-known union bound method 
of H. It is also shown that this upper bound can be computed 
with a closed form expression in the case of real constellations 
and uncorrelated fading conditions. However for complex 
constellations and over correlated fading channels, it needs 
to be computed numerically using the Monte Carlo methods. 
One alternative to numerical averaging methods is to employ 
the Chernoff or other exponential upper bounds as done in 
Bl , |6), but these bounds are not tight especially in the case 
of correlated channels and higher spectral efficiencies. Notice 
that the performance analysis of SSK poses a relatively simpler 
problem compared to SM since only the antenna indices are 
transmitted as information. That is why, accurate closed form 
upper bounds are derived for Rayleigh fading channels in 0, 
for correlated Rician channels in [7| and for single receiver 
Nakagami fading channels in (8). Because of its inherent 
difficulty in the general MIMO setup, the SM performance 
analysis problem is studied with only one receive antenna 
in 151 for uncorrelated and in iflOl for correlated Rayleigh 
fading channels. Among these, is particularly important 
as the authors make a detailed analysis of the error metric 
and show that by employing a Gaussianity approximation it 
is possible to derive tight upper bounds for certain digital 
modulation schemes such as phase-shift keying (PSK) and 
multilevel quadrature amplitude modulation (QAM). However, 
the analytical results of [9] and [ 1 1 are not extendable to 
general SM schemes with multiple receive antennas and/or 
with other modulation schemes such as rectangular QAM. 

To the best of our knowledge, performance analysis of SM 
with an arbitrary number of transmit/receive antennas and 
general family of linear modulation alphabets is still an open 
problem. In the error analysis of conventional communication 
systems, the envelope of channel fading effects (in scalar or 
vector form) appears in the error performance metric. The 
distribution of this envelope or its square can be obtained 
for most known fading conditions. Then using probability 
density function (PDF), more commonly moment generating 
function (MGF) (Til. C2 or characteristic function (CHF) 



Ifl3l approaches it is possible to compute the average error rate 
for a wide class of digital modulation methods. As pointed out 
El, 0, [TO), the difficulty in SM arises from the fact that the 
error metric consists of the envelope of not a single fading 
variable. Instead it contains a complex weighted mixture of 
two complex random vectors formed by the channel fading 
coefficients where the mixture weights are the constellation 
symbols. Even in the uncorrelated Rayleigh fading case, if the 
constellation points are not real, the real and imaginary parts 
of the resulting vector are not statistically independent. This 
makes it challenging to obtain a PDF or MGF (or CHF) for 
the envelope and therefore to derive a closed form expression 
for the average pairwise error probability (APEP). 

In this paper, we propose an upper bounding technique 
for SM that is applicable to an arbitrary number of trans- 
mit/receive antennas and general modulation alphabets for 
uncorrelated/correlated Rayleigh and Rician fading channels. 
We extend the insight of [2| that is limited to SM with real con- 
stellations to the case of complex constellations. Specifically, 
we show that for a number of fading scenarios, the random 
variable/vector that appears in the SM error metric is a "proper 
complex random variable/vector" defined in |[T4l and can be 
described with a joint multivariate Gaussian PDF. Then using 
its multivariate distribution instead of that of its envelope and 
employing the approach in [15], we present a framework to 
compute the upper bound for a class of fading channel models. 
We assume a general parametric correlated Rician channel 
model that can also be used to describe uncorrelated Rician 
and Rayleigh fading conditions. The general framework can 
be utilized to derive error upper bound for a large class of 
channel fading conditions. The framework provides an exact 
bound generally via simple numerical integration and in some 
cases results in closed-form bounds. 

The rest of this paper is organized as follows. In Section II, 
the system and generalized channel fading model is presented. 
Then a brief review of proper complex random variables is 
given and the performance bounds are derived for MIMO SM 
systems in Section III. This paper is ended with conclusive 
remarks and discussions in Section IV. 

II. System Model 

A. Signal Model 

We consider the conventional SM system model with Nt 
transmit and N r receive antennas employing optimal detec- 
tion in J2). We also assume that the number of transmit 
antennas is an integer power of 2, i.e., Nt = 2™, and the 
transmitter employs M-ary digital modulation to an m-bit 
message where M = 2 m and the modulation symbol set 
is X = { Xi, ... Xk, ••• X M }. Notice that the 
spectral efficiency of the system is R = n + rn = \og 2 (N t M). 
At each transmission instant, each set of n + to of bits is 
split into groups of n and to bits, and the former is used to 
select one of the Nt antennas and the latter to be mapped onto 
one of M possible complex constellation points determined by 
the particular digital modulation method. The method used 
in the bit-to-antenna index and bit-to-symbol mappings is 



immaterial to the central discussion, so throughout the paper 
we assume uniform mapping. The transmitted signal vector is 
x = [ x\, ... xi, ... xjsr t 1 where all but one entry 
is zero because only one antenna is active for transmission. 
Notice that if the £-th antenna is selected, then all entries other 
than X£ is zero and xe £ X. In other words, the position of 
the non-zero element denotes the antenna index and its value 
indicates the transmitted symbol. Similar to 0, we assume 
a power constraint of unity, i.e., E x [x^x = l] . The received 
signal model is expressed as 



VpHx 



(1) 



where p = E s /Nq is the average signal to noise ratio (SNR) 
observed at each receiver branch, y and v are the N r X 1 
received signal and channel noise vectors, respectively, and H 
is the N r x Nt dimensional channel matrix. The elements of 
v are modelled as independent identically distributed (i.i.d.) 
complex Gaussian variables with zero mean and unit variance, 
i.e., v k ~ C7V(0, 1) for k = 1, . . . , N r . 

B. Fading MIMO Channel Model 

We assume a slow fading MIMO channel model with the 
sum of an average (or fixed, possibly line-of-sight) component 
and a variable (or random) component. Accordingly, N r x Nt 
dimensional channel matrix H is described as 
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(2) 



where H and H are the fixed and variable components, respec- 
tively. The square-root terms are the normalization weights 
with K being the Rician factor and reflecting the ratio of the 
fixed and variable channel components. The fixed part H is 
modelled as an all-one matrix, i.e., [H] Uj „ = h u ,v — 1, for 
u = 1, . . . , N r and v — 1, . . . , Nt. The variable component 
of the channel matrix, H, consists of -possibly correlated- 
complex Gaussian variables. Given [H] Uj „ = h u ,v, f° r all 
channel coefficient pairs (h u>V) h& t %) (u,ii = l,...,N r and 
v, v = 1, . . . , N t ) we assume that 
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that is to say, the auto-correlations of the real and imaginary 
parts are the same and there is no correlation between real and 
imaginary parts. With this condition, the correlated channel 
matrix can be described by the well-known Kronecker corre- 
lation model in which H is expressed as 



H = sy 2 Hsf /2 



(3) 



where S t and S r are the real valued and Hermitian symmetric 
transmit and receive correlation matrices, respectively, with 
the elements defined as [J7 t ] Uj u — a l u A for u,u= 1, . . . , Nt, 
[£ r ] v $ = o r v ~ for v, v = l,...,N r . H is the independent 
Rayleigh fading channel matrix with zero mean and unit 
variance. In other words, the elements of H are described as 



independent identically distributed complex Gaussian random 
variables, i.e., [H] u .„ = h u . v ~ CAf(0, 1) for u = 1, . . . , N r 
and v — 1, . . . ,N t . Notice that, with this channel model the 
multivariate conditional PDF of y can be written as: 

/ Y (y |x,H)= ^^^(s) exp ("(y- m Y) t S z ~ 1 (y-m Y )) 

(4) 

where m Y = E [y] = J y^-Hx and S = 

are the mean 



: E[y] = ^Hx and 
E[(y-m Y )(y-m Y )t] = ® S. 

vector and the covariance matrix, respectively (The operator 
® denotes the Kronecker product). Notice that the Kronecker 
correlation model separating the transmit and receive correla- 
tion effects has been used extensively within the context of 
SM/SSK literature JH-Q. The correlation matrices can be 
formed according to a number of models as detailed in the 
discussion of the simulation parameters and results. 

Combining (|2} and (01, the general channel model can be 
rewritten as 



which, via the appropriate selection of parameters, makes it 
possible to characterize uncorrelated/correlated Rayleigh and 
Rician fading scenarios. 

C. Optimal Detector and Average Probability of Error 

Given the conditional PDF in ©, the optimum antenna and 
symbol index pair (u]yiL> w Ml) m tne maximum likelihood 
(ML) sense is expressed as 

( W ML' U ML) = ar g max /Y(y | x,H) = argmini?(y,h u X t ,) 

u,v u : v 

where h u is the w-th column of H for u = 1, . . . , N t and 
X v is the v-th element of the modulation alphabet X for 
v = 1, . . . , M. D(y, h u X v ) is the -modified- distance metric 
between y and h u X v defined as 

D(y, h u X v ) = J-p || h u X v || 2 -2Re{yth u Xj. (6) 

As shown in [2|, the average bit error performance of the 
SM system with optimum detection can be computed by the 
well-known upper bounding technique in [4|. Accordingly, the 
average bit error probability (ABEP) is upper bounded as 

N t N t M M 

Pb < 



In ([U, the N r x 1 vector z is defined as J2): 

z = J~^(h u X v - huXv). 

Notice that the fc-th element of z, k = 1, . . 
defined as 



(9) 



. , N r can be 
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(10) 

z k = c[(X R -X R ) + i(X^-Xi)] (11) 
Zk — c {h R u X R — hf, u X* — h R ^X R + h{^Xf^j 

+ * c [hk,uX v + h{ u X R — hf.^Xjj — h{ fi Xjfj(12) 



1 

N t M 



iLLLL ^/» r ^ Ps(u,u,v,v) (7) 



log 2 (N t M) 



u—l u—1 v—1 v — 1 

where N(u, u, v, v) is the number of bits in error between the 
respective channel and symbol pairs, (h u ,X v ) and (h^, X v ). 
The denominator term \og 2 (N t M) in expression in rep- 
resents the total number of antenna and symbol bits and 
division with this term indicates the summation weight for 
the corresponding pairwise error probability (PEP). For large 
spectral efficiencies (bits per SM symbols) the average of the 



where the coefficient c = y 2(if+i) ' an< ^ Zk an< ^ * k denote me 
fixed and variable components of Zk, respectively. 

Notice that in 0, it is shown that so long as the symbol 
constellation X is real, the real and imaginary parts of the 
expression in ( fT2l are statistically independent of each other 
in the case of uncorrelated Rayleigh fading channels (K = 0, 
H = H). Then for any symbol pair (X v , X$), the term |j z || 2 
becomes a Chi-squared random variable with 2N r degrees 
of freedom and therefore the average of the Q-function in 
dT) can be computed in closed form. On the other hand, as 
also indicated in [2|, when (X V ,X V ) come from complex 
constellations the distribution of the envelope of random 
variable Zk (or z& in general) cannot be easily obtained, due 
to statistical dependency of the real and imaginary parts of Sk- 
in this case, the Monte Carlo averaging methods are adopted 
in both 12) and the subsequent works (9), (TO) for both SM 
and SSK error analysis. 

In the next section, we show that accurate statistical descrip- 
tion of z is possible in terms of a joint multivariate PDF so 
long as it is a proper complex Gaussian random vector. Then 
using this distribution of z we provide a general framework in 
which the APEP in <(8J and the ABEP in (O can be computed. 

III. Performance Analysis of Spatial Modulation 

A. Proper Complex Random Vectors 

Let Y = Y R +iY I be a N x 1 complex random vector with 
the mean m Y = E [Y] . Suppose its real covariance matrices 
Sy, ^y 1 an( l 27 Y an d com pl ex covariance and pseudo- 
covariance matrices Sy and Sy are respectively given by the 
following equations as 



ratio 



N (u,u,v ,£>) 
log 2 (JV t M) 



approaches to the factor of |, which is the one 
used in the derivations [9], [10|. In (0, P s (u,u,v,v) is the 
average pairwise symbol error probability (APEP) as defined 
in equation ([SJ. 
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P s (u, U, V, ' 



E x ,h \Pr{D(y,h u X v ) > Z?(y,h^X 5 ) | 



(8) 



Definition [14|: Y is said to be proper complex random 
vector if its pseudo-covariance matrix Sy = 0, or equiva- 
lently, if E$ = Z{- and S^ 1 = S^ R . 

A direct consequence of this condition is that 



Sy = 2Si 



i2S 



in 



(13) 



In the special case of Y being a scalar, denoted by Y, 
E Y R — S RI . If Y is proper, which means E [(Y — m Y ) 2 ] = 
0, then Sy — —S RI together with the condition above 
implies that U IR = 0, i.e., Y R and Y 1 are uncorrelated. 
This along with d!31 l. implies that for a complex scalar Y 
to be proper the real and imaginary components must have 
the same variance and be uncorrelated. 

Finally if Y is a proper complex Gaussian vector, Y R and 
Y 7 are jointly Gaussian and the joint multivariate PDF of Y 
is written as 

= ^d^) e " Cy_mY)t ^ (y_mY) - (U) 

Given the definitions of proper random vectors above and 
z in (|9|>-([T2||. it is easy to show that z is a proper Gaussian 
random vector and derive the parameters of its PDF for both 
N r = 1 and N r > 1. 

When N r = 1, notice that the real and imaginary parts of 
z are Gaussian random variables (real weighted sum of real 
Gaussian variables). Furthermore, 



m z 
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In other words, the real and imaginary parts of z are uncor- 
related and the pseudo-variance of z vanishes. This makes 
the variance real and therefore z a proper complex Gaussian 
random variable by definition. 

When N r > 1, notice from above that each element of the 
vector z consists of a proper complex Gaussian random vari- 
able with uncorrelated real and imaginary parts and identical 
-real- variances. Given the receiver correlation matrix 27 r , and 
the definition of z we can show that: 
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c 2 [| X v | 2 + | Xt, | 2 -2Re{<7* ie X„X e *}] S T 



Z RI = 



c 2 [| X v | 2 + | Xd I 2 -2Re{al & X v X;}] 27 r (16) 



This implies that the vector z is a proper complex Gaussian 
vector with the conditional joint PDF given as 



/z(z | h u ,hu,X v ,Xi,) = 



-(z-m z ) t i:- 1 (z-mz) 

Tr^det^z) 



(17) 



where det(-) represents the matrix determinant operation and 
mz and Sz are as given in (l5[ and ( [Tol l, respectively. Notice 
that proper complex random variables and vectors whose 
properties are presented in [14] do not require the real and 
imaginary components to be statistically independent. So long 
as the pseudo-covariance matrix vanishes, a proper complex 
random variable (vector) has a distribution. 

B. Upper Bound for the Error Performance of SM 

The distribution obtained in ( TTTb can be used to evaluate 
the APEP in ([8]) by following an approach similar to that of 
IfTSI , such that 

P s (u,u,v,v) = 

^Q(Vl|z|| 2 )/z(z | h„,h fi) X^,X«)dz.(18) 

Then by using the properties of z and Craig's alternative 
definition of the Q(-) function, (fT8l is rewritten as in the 
integral formulations given the top of this page where the 
parameter change z z/ v2 results in 
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Notice as in the second line of ( lfl4l . eq. (6)), by the Hermitian 
symmetry property of Sz and the fact that the joint PDF 
of a proper complex Gaussian vector integrates to unity, this 
integral then simplifies to the expression below 



P„(u,u,v,v) 



1 



[-Ez+sin 2 01]- 



det 



+ 1 ) 



-de. 



(19) 



Plugging this result in (|7]) gives the ABEP expression for SM 
employing any modulation alphabet. 

The APEP integral in ( fT~9b simplifies to some closed form 
expressions in certain cases as we show in the sequel. Other- 
wise, especially in the general Rician fading case, i.e. when 
mz 7^ 0, no further simplification is possible but since < TT~9b 
is a single definite integral, it can be conveniently evaluated 
numerically for any given mean vector and covariance matrix 
without the need for any Monte Carlo averaging. 
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C. Special Case: Average Error Performance for SSK 

Notice from the definition of z in (0 that under the special 
conditions that h u = ha (same transmit antenna), X v = X{, 
(same transmitted symbol), z remains to be a proper complex 
random vector. Of these, the latter is particularly important 
as it corresponds to the transmission of information bits with 
single symbol constellation, or equivalently, with SSK. Notice 
that in this case z becomes 

z = ^(h u X v - h fi X„) =c(K- h«) X v (20) 

with the corresponding first and second order statistics 
m z = 

*z = 5 ^ iy (2-2Re[< 4 J)^. 

The condition in (l20b implies that so long as the fixed 
components of the MIMO Rician channel are the same, their 
contributions cancel out in the expression of z and the error 
performances over correlated Rayleigh and Ricean channels 
differ only with the scaling factor of \j\/K + 1. Accordingly, 
the ABEP for SSK can be written as 

N t N t 



pSSK < 



N, 



J_ sr v NM p ssk (u u) 



where 



Pf SK (u,u) = - 
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det 
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with Si 



SSK 



'z ^ 2(rffTJ ( X - Re [<fil) As a result the 
error upper bound derivations of SSK can be computed as 
that of SM with M-PSK with M = 1. 

IV. Simulation Results 

In this section, simulation results are presented to validate 
the proposed error performance calculation method for the 
SM/SSK systems under various channel fading conditions. For 
all simulations, a MIMO system with 4 transmit and 4 receive 
antennas is considered. We simulate transmission schemes 
with spectral efficiencies of R = 3 b/s/Hz to R = 7 b/s/Hz. In 
a Nt = 4 system this corresponds to employing modulation 
alphabets of 2 up to 32 symbols. For small R, BSPK and 
QPSK constellations are used while for large R values they 
are replaced by rectangular QAM constellations with 8, 16 
and 32 points. Notice that rectangular QAM constellations are 
chosen in order to show that the approach presented in this 
paper is not dependent on a particular constellation geometry. 
In all of the results presented, the curves corresponding to 
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Figure 1: SM error performance with Nt = 4 and N r = 4 
fori? = 3-7 b/s/Hz. Uncorrelated Rician fading with K=5. 



simulations are plotted with solid lines whereas those from 
analytical upper bound derivations are represented with dotted 
lines with identical markers. For reasons of clarity the marker 
of only one curve (simulation) is included in the legends. 

The transmit and receive channel correlation matrices in 
the Kronecker model are formed according to two common 
and frequently used models. In the first, presented in fl6l . 
the correlation matrices are computed based on a clustered 
channel model using the power azimuth spectrum distribution 
and the array geometry. The exponential correlation model of 
ifTTll is also considered where the correlation matrix entries are 
formed as [27] = a u _ v = j\ u ~ v \ where f3 is a fixed (real or 
complex) correlation coefficient between adjacent antennas. In 
simulations, we consider the cases where i) jt = 7r = 0.5, ii) 
j t = 7 r = 0.8, and hi) j t = 0.9, 7 r = 0.1. Notice that the 
first and second cases are used to model moderate and strong 
correlation on both sides, the third to describe heavy transmit 
correlation with almost no correlation at the receiver side. 

The theoretical error upper bound is compared with bit error 
performances obtained through numerical simulations in Fig. 
03 -Fig- El Fig- 03 shows performance of SM for uncorrelated 
Rician fading channels with a Rician factor of K = 5 whereas 
the performance of SM under correlated Rician channels is 
shown in Fig. [2] Finally, Fig. [3]presents the SSK performance 
comparisons for Nt — 4 (R — 2 b/s/Hz spectral efficiency) 
under correlated Rician fading. Notice that in all SM simu- 





i i i ; i : ; ; 




■ .:■ :■ ! ; d- : 
— *— 4 h/s/Hz: OPSK 
— — 5 bls/Hr. B-QAM - 


r 
























; 












' ' ' 5 « 




: 






:: :: 5:: :: I ^ !i :: 


% W 




: 






X 


I 























6 


12 15 18 21 i 


:;:::':-;x::- 

4 27 30 33 36 39 42 



Figure 2: SM error performance with N t = A and iV r = 4 for 
i? = 3-7 b/s/Hz. Correlated Rician fading, K=5. Exponential 
correlation model of ifTTll with parameters j t — 7r = 0.8. 




Figure 3: SSK error performance with Nt = 4 and N r = 4 
for correlated Rician fading, K=5. 

lations the results are depicted for the exponential correlation 
model [17 1 where the clustered correlation model is omitted 
for space limitations. Both correlations models are depicted in 
the SSK simulations. All curves are drawn down to at least 
1 x 10~ 6 bit error rate (BER) levels and simulations indicate 
the upper bounds obtained with analytical framework are tight 
in all channel conditions considered. 

V. Conclusion 

In this paper, we present a performance analysis framework 
for SM schemes over correlated Rayleigh and Rician fading 



channels characterized by Hermitian symmetric real corre- 
lation matrices. Our framework is based on the observation 
that the complex vector used in the error function is proper 
complex Gaussian and has a multivariate joint distribution. We 
show that using this distribution it is possible to obtain exact 
bounds on the performance of SM schemes and a number of 
special cases. One particular case is SSK, where the upper 
bound can be computed readily by considering it as SM with 
a constellation size of 1. This upper bounding method for SSK 
is consistent with the physical realization of the transmission. 
Simulation results indicate a tight match with the theoretical 
derivations. 
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